Evaluation of production performance from a hydraulically fractured well

ABSTRACT

An analytical solution is obtained for a pseudo-steady state production from a vertically fractured well with finite or infinite fracture conductivity. The analytical solution may be used to compute a pseudo-steady state constant for the reservoir. Subsequently, performance parameters relating to the reservoir may be derived from the pseudo-steady state constant. For example, parameters such as production decline rate, total hydrocarbon reserves, and economically recoverable reserves for the reservoir may be computed. The disclosed analytical solution, instead of a conventional numerical simulation, can significantly speed up analysis and improve the accuracy of the calculation of these parameters.

CROSS-REFERENCE TO RELATED PATENT APPLICATIONS

This application claims the benefit of U.S. Provisional Patent Application No. 62/268,958 to Kangping Chen, filed on Dec. 17, 2015, and entitled “Evaluation Of Production Performance From A Hydraulically Fractured Well,” which is hereby incorporated by reference in its entirety.

FIELD OF THE DISCLOSURE

The instant disclosure relates to extraction of underground resources. More specifically, this disclosure relates to determining performance factors relating to the extraction of underground resources from a particular well.

BACKGROUND

Production of hydrocarbon from a well is normally conducted with a constant production rate over long periods, although the rate can be changed during the productive life of the well due to maintenance and other technical requirements. FIG. 1 is a graph illustrating reservoir pressure change with time for a well producing at constant rate in a closed reservoir. P_(i) is a reservoir initial pressure; P_(w) is a wellbore pressure; and P_(cri) is the lowest permissible wellbore pressure (critical pressure). The time sequence of the graph of FIG. 1 is t₁<t₂<t₃<t₄<t₅<t₆ . . . . At the start of a production, reservoir pressure initially depletes in the immediate neighborhood of the wellbore, and this pressure drawdown spreads outward diffusively towards the reservoir outer boundary (as shown in FIG. 1). For a closed (sealed) reservoir, the no-flow reservoir boundary starts to affect the pressure when the spreading pressure depletion front approaches the boundary. When the boundary effect has been fully reflected in the pressure field, the spatial distribution of the pressure no longer changes with time and the fluid flow reaches the so-called pseudo-steady state (lines 102A in FIG. 1). The flow prior to the pseudo-steady state flow is called the transient flow (lines 102B in FIG. 1), the duration of which depends on how fast the pressure drawdown diffuses in the reservoir, which in turn is determined by the reservoir and fluid properties, namely permeability, porosity, viscosity and compressibility. For conventional reservoirs where the permeability is greater than 0.1 mD (mini-Darcy), the transient flow period usually lasts from days to months; while for unconventional reservoirs which have permeabilities less than 0.1 mD, the period can last from years to even tens of years. For closed reservoirs, the pseudo-steady state flow is a dominant, long-duration and most productive flow regime, especially for conventional reservoirs. During the pseudo-steady state flow period, the wellbore bottom-hole flowing pressure (BHFP) decreases linearly in time in order to maintain the constant production rate. However, once the bottom-hole flowing pressure has declined to the lowest permissible value, which is often determined by the surface equipment limitations, a constant rate production can no longer be continued, and a constant pressure production must follow. The production rate for this constant pressure production period declines in time, eventually approaching zero as the reservoir pressure approaches the lowest permissible wellbore pressure (lines 102C in FIG. 1).

Pseudo-steady state flow is a dominant flow regime during constant rate production from a finite, closed reservoir. For a vertically fractured-well in a finite reservoir approximated as having an slightly elliptical shape, conventional solutions exist for analytically determining the flow for the case of infinite fracture conductivity. For finite fracture conductivity, conventional computational techniques to achieve a pseudo-steady state solution involve running numerical simulations over long times of hours, days, or longer.

Pseudo-steady state flow is a dominant flow regime during constant rate production from a closed reservoir: after the effects of the no-flow condition on the reservoir outer boundary have been fully reflected in the flow field and the transients associated with the flow startup have decayed to be negligible, the flow in the reservoir reaches a state in which the spatial distribution of the pressure no longer changes with time. Pseudo-steady state flow is thus a boundary-dominated flow. One definition for pseudo-steady state is the condition in a finite, closed reservoir when producing at a constant rate that “every point within the reservoir will eventually experience a constant rate of pressure decline.” This constant rate of pressure decline is the result of mass conservation for constant rate production from a closed reservoir. This condition is sometimes referred to as pseudo-steady, quasi-steady, semi-steady, or even steady state. The term pseudo-steady is used here in reference to this particular flow regime.

Pseudo-steady state (PSS) can be a prolonged period of constant rate production from a closed reservoir. During this period, the reservoir pressure declines linearly with time, the rate of which is determined by the specified production rate and the drainage area. The pseudo-steady state solution provides the reservoir pressure distribution as well as the productivity index for this important flow period. Once the bottom hole flowing pressure has declined to the lowest permissible value, however, a constant rate production can no longer be continued, and a constant pressure production must follow. The production rate for this latter constant pressure production period declines in time. Production rate decline analysis for this period plays an important role for estimating the hydrocarbon reserves in place and for assessing the economically recoverable amount of fluid from a reservoir. Because pseudo-steady state is the flow regime immediate preceding the production rate decline period, the pseudo-steady state solution has been conventionally used in the production rate decline analysis for unfractured wells and for fractured wells. In these analyses, the pseudo-steady state dimensionless pressure drawdown at the wellbore is expressed as Δp _(wD,PSS)=2πt _(DA) +b _(D,PSS),  (1) where t_(DA) is the drainage area based dimensionless time, and b_(D,PSS) is the so-called pseudo-steady state constant which depends on the reservoir model as well as the well/reservoir configuration. This pseudo-steady state constant b_(D,PSS) is used to define the appropriate dimensionless decline rate and time in many of the currently used production decline rate analysis models. Furthermore, the pseudo-steady state constant is the reciprocal of the dimensionless productivity index J_(D,PSS) for the pseudo-steady state, J_(D,PSS)=1/b_(D,PSS), which measures the productivity of the well for this flow period. J_(D,PSS) is also important for production optimization for a fractured well. For unfractured wells, the pseudo-steady state constant b_(D,PSS) can be obtained analytically for reservoirs of very simple shapes. These exact analytical solutions have been modified by shape factors and used as approximate analytical solutions for other reservoir geometries. For hydraulically fractured wells, however, exact analytical solution for b_(D,PSS) is not available. For a vertically fractured well with infinite fracture conductivity, an exact analytical solution for the pseudo-steady state flow in a reservoir bounded by an elliptical boundary is known, which leads to an analytical expression for the pseudo-steady state constant b_(D,PSS). For the more practical case of finite fracture conductivity, however, no exact analytical solution in the physical variable space has been reported in the literature for pseudo-steady state flow. For finite fracture conductivity, one conventional numerical procedure is to extract b_(D,PSS) by subtracting 2πt_(DA) from the long-time numerical solution for constant rate production from a fractured well in an elliptical reservoir. This procedure is quite time consuming; and curve-fitting has been used to obtain an empirical relation between b_(D,PSS) and the reservoir geometric parameter and the fracture conductivity.

SUMMARY

An analytical solution for pseudo-steady state flow for a vertically fractured well with finite fracture conductivity in a closed reservoir modeled as having a nearly circular, slightly elliptical shape is described in embodiments of the present invention. This analytical solution provides a solution to a problem with no previous known analytical solution. Furthermore, the analytical solution can be used in computer simulations to improve production performance of a hydraulically fractured well, provide prospectors with improved information for deciding on production wells, and improve production from those wells selected for production. The analytical solution allows computer modeling to be performed accurately and timely. Conventional techniques described above failed to provide an analytical solution for pseudo-steady state flow for vertically fractured wells, and those conventional techniques consumed significant amounts of computer processing time.

The analytical solution can be expressed in terms of elementary functions and provides a simple expression for the pseudo-steady state constant and the dimensionless productivity index. This analytical solution may be executed on a computer system to quickly generate performance parameters or other characteristics of the vertically fractured well. This solution eliminates the need of performing time-consuming numerical simulation for obtaining pseudo-steady state solution for fractured wells in a near circular reservoir and it may be used to generate approximate solutions for reservoirs of other geometrical shapes. For example, in comparison to the hours or days required of a computer to generate solutions according to the conventional techniques described above, a computer may generate solutions in accordance with described embodiments of the invention in a matter of seconds or minutes.

Described embodiments may yield a simple, exact expression for the pseudo-steady state constant b_(D,PSS), which can be used for various applications including production rate decline analysis and fracture design for optimized production. The solution can also be used as a benchmark to measure the accuracy of numerical simulations. With suitable shape factors, the analytical solution may be used to obtain approximate expressions for the pseudo-steady state constant b_(D,PSS) for fractured wells in reservoirs of other geometrical shapes.

According to one embodiment, a method may include receiving a plurality of shape factors corresponding to a geometrical shape of a hydraulically fractured well reservoir; determining a pseudo-steady state constant for the reservoir based, at least in part, on the plurality of shape factors; and/or determining a performance parameter of the reservoir when operated in a pseudo-steady state with a finite fracture conductivity based on the determined pseudo-steady state constant.

The foregoing has outlined rather broadly the features and technical advantages of the present invention in order that the detailed description of the invention that follows may be better understood. Additional features and advantages of the invention will be described hereinafter that form the subject of the claims of the invention. It should be appreciated by those skilled in the art that the conception and specific embodiment disclosed may be readily utilized as a basis for modifying or designing other structures for carrying out the same purposes of the present invention. It should also be realized by those skilled in the art that such equivalent constructions do not depart from the spirit and scope of the invention as set forth in the appended claims. The novel features that are believed to be characteristic of the invention, both as to its organization and method of operation, together with further objects and advantages will be better understood from the following description when considered in connection with the accompanying figures. It is to be expressly understood, however, that each of the figures is provided for the purpose of illustration and description only and is not intended as a definition of the limits of the present invention.

BRIEF DESCRIPTION OF THE DRAWINGS

For a more complete understanding of the disclosed system and methods, reference is now made to the following descriptions taken in conjunction with the accompanying drawings.

FIG. 1 are graphs illustrating reservoir pressure change with time for a well producing at constant rate in a closed reservoir according to the prior art.

FIG. 2 is a top view of a vertical well intersected by a thin elliptical fracture according to some embodiments of the disclosure.

FIG. 3 is a flow chart illustrating an example method for computing an analytical solution for pseudo-steady state flow for a vertically fractured well with finite fracture conductivity in a closed reservoir according to some embodiments of the disclosure.

FIG. 4 are graphs of a pseudo-steady state constant computation as a function of ξ_(e) calculated according to some embodiments of the disclosure.

DETAILED DESCRIPTION

Consider fluid production from a fully-penetrated, vertically-fractured well from an initially quiescent state as shown in FIG. 2. FIG. 2 is a top view of a vertical well intersected by a thin elliptical fracture according to some embodiments of the disclosure. The drawing is for illustration purpose only and it does not reflect the actual scales. In some embodiments, the fracture may be very thin and long, L>>w_(f), and the fracture surface, ξ=ξ₁≈0. The following commonly used assumptions are made: the reservoir fluid is a single phase fluid residing in a homogeneous medium with its motion governed by the Darcy's law in both the reservoir and the fracture; the fluid and the reservoir are weakly compressible, characterized by a single lumped total compressibility constant c_(t); the effects of wellbore storage and skin are negligible; and the hydraulic fracture is supported by propants and it is incompressible. The hydraulic fracture is modeled as a thin, long ellipse, intersecting the wellbore with a fracture width w_(f), which is much smaller than the wellbore diameter 2r_(w). The Cartesian coordinates (x,y) and the elliptic coordinates (ξ,η) are related by x=L cos hξcos η, y=L sin hξ sin η, with L being the focal distance which is essentially the fracture half-length; and L>>w_(f). The surface of the narrow elliptical shape fracture is represented by the ellipse ξ=ξ₁ in the elliptic coordinates, and ξ₁ is a small number. Subscript “f” is used for reservoir and fracture quantities, respectively. The permeabilities in the reservoir and the hydraulic fracture are κ,κ_(f), respectively, with κ_(f)>>κ. For a successful hydraulic fracturing job, the fluid production is nearly entirely from the fracture, and the contribution from the wellbore to the production is negligible. The reservoir has a finite extent and its outer boundary is an ellipse ξ=ξ_(e), confocal with the limiting ellipse ξ=ξ₁ used to represent the fracture. For mathematical simplicity, the finite drainage area is assumed to have an elliptical shape, which is a good geometrical approximation to a large circular drainage area. A large circular reservoir with radius R can be well approximated as an elliptical reservoir with ξ_(e)=ln(2R/L).

For the convenience of discussing the physical aspects of the analytical solution, we formulate the problem in terms of pressure instead of pressure drawdown as in most of petroleum engineering literature. Pressure drawdown will be denoted as Δp throughout the paper. We choose the reservoir initial pressure p_(i,d) and the pressure diffusion time scale as the characteristic pressure and characteristic time, respectively, for non-dimensionalization: p_(c)=p_(i,d),t_(c)=μφc_(t)L²/κ, where μ,ϕ are the fluid viscosity and reservoir porosity, respectively. The dimensionless reservoir pressure satisfies a diffusion equation, which in an elliptical coordinates (ξ,η) becomes

$\begin{matrix} {{{\frac{\partial^{2}p_{D}}{\partial\xi^{2}} + \frac{\partial^{2}p_{D}}{\partial\eta^{2}}} = {\frac{{\cosh\; 2\;\xi} - {\cos\; 2\;\eta}}{2}\frac{\partial p_{D}}{\partial t_{DL}}}},} & (2) \end{matrix}$ where the dimensionless time and the dimensionless pressure are defined by t _(DL) =κt/(μφc _(t) L ²),  (3) p _(D) =p _(d) /p _(i,d).  (4) Initially the reservoir fluid is at rest. Symmetry condition applies on the x-axis and the y-axis; and the no-flow condition is imposed on the reservoir outer boundary,

$\begin{matrix} {\xi = {{\xi_{e}\text{:}\frac{\partial p_{D}}{\partial\xi}} = 0.}} & (5) \end{matrix}$ When the fracture compressibility is neglected, the dimensionless pressure in the fracture, defined as p _(fD)(η,t _(DL))=p _(f,d)(η,t _(DL))/p _(i,d),  (6) satisfies the equation

$\begin{matrix} {{{{\frac{\partial^{2}{p_{fD}\left( {\eta,t_{DL}} \right)}}{\partial\eta^{2}} + {\frac{2}{F_{E}}\frac{\partial p_{D}}{\partial\xi}}}}_{\xi\; = \;\xi_{1}} = 0},} & (7) \end{matrix}$ where the dimensionless elliptical fracture conductivity

$\begin{matrix} {F_{E} = {\frac{\kappa_{f}w_{f}}{\kappa\; L}.}} & (8) \end{matrix}$ This elliptical fracture conductivity F_(E) is different from the rectangular fracture conductivity commonly denoted as C_(fD). For an elliptical fracture, the width of the fracture is not a constant; F_(E) and C_(fD) only match with each other at the well. One way to relate F_(E) and C_(fD) is to assume that the elliptical fracture and the rectangular fracture have the same volume, which leads to C_(fD)=πF_(E)/2. Symmetry condition on the x-axis is imposed and a constant production rate at the wellbore is specified. The Pseudo-Steady State Pressure Distribution in a Closed Elliptical Reservoir

A pseudo-steady-state (PSS) solution is the long-time asymptotic solution under constant production rate condition from a closed reservoir; and it has the property that

$\begin{matrix} {{\frac{\partial p_{D,{PSS}}}{\partial t_{DL}} = {{{const}.} = {{- C} < 0}}},} & (9) \end{matrix}$ where C is a dimensionless positive constant, C>0. Property (9) holds for any point in the reservoir. Thus, the reservoir pressure possesses the form p _(D,PSS)(ξ,η,t _(DL))={tilde over (p)}(ξ,η)−Ct _(DL),  (10) and eqn. (2) becomes an eqn. for the shape function {tilde over (p)}(ξ,η):

$\begin{matrix} {{\frac{\partial^{2}\overset{\sim}{p}}{\partial\xi^{2}} + \frac{\partial^{2}\overset{\sim}{p}}{\partial\eta^{2}}} = {{- C}{\frac{{\cosh\; 2\;\xi} - {\cos\; 2\;\eta}}{2}.}}} & (11) \end{matrix}$ The solution to the inhomogeneous eqn. (11) can be written as {tilde over (p)}(ξ,η)=p _(c)(ξ,η)+p _(p)(ξ,η),  (12) where p_(c)(ξ,η) satisfies the homogeneous eqn.

$\begin{matrix} {{{\frac{\partial^{2}p_{c}}{\partial\xi^{2}} + \frac{\partial^{2}p_{c}}{\partial\eta^{2}}} = 0},} & (13) \end{matrix}$ and p_(p)(ξ,η) is a particular solution of the inhomogeneous eqn.

$\begin{matrix} {{\frac{\partial^{2}p_{p}}{\partial\xi^{2}} + \frac{\partial^{2}p_{p}}{\partial\eta^{2}}} = {{- C}{\frac{{\cosh\; 2\;\xi} - {\cos\; 2\;\eta}}{2}.}}} & (14) \end{matrix}$ One solution to equation (14) is

$\begin{matrix} {{p_{p}\left( {\xi,\eta} \right)} = {{- \frac{C}{8}}{\left( {{\cosh\; 2\xi} + {\cos\; 2\eta}} \right).}}} & (15) \end{matrix}$ A possible solution to the homogeneous equation (13) has the form:

$\begin{matrix} {{{p_{C}\left( {\xi,\eta} \right)} = {B_{0} + {A_{0}\xi} - {\sum\limits_{n = 1}^{\infty}{A_{n}\cos\; 2n\;{\eta cosh2}\;{n\left( {\xi_{e} - \xi} \right)}}}}},} & (16) \end{matrix}$ where A_(i)(i=0, 1, 2 . . . ) are constants, and the symmetry conditions on the x-axis and y-axis (η=0,π/2) have already been satisfied. The infinite series enters the fracture eqn. because of its non-zero flux density on the fracture surface. For the case of finite fracture conductivity, this infinite series is needed to match the non-constant fracture pressure inside the fracture.

Thus, the reservoir pressure for pseudo-steady state is given by

$\begin{matrix} {{p_{D,{PSS}}\left( {\xi,\eta,t_{DL}} \right)} = {B_{0} + {A_{0}\xi} - {\sum\limits_{n = 1}^{\infty}{A_{n}\cos\; 2n\;{\eta cosh2}\;{n\left( {\xi_{e} - \xi} \right)}}} - {\frac{C}{8}\left( {{\cosh\; 2\xi} + {\cos\; 2\eta}} \right)} - {{Ct}_{DL}.}}} & (17) \end{matrix}$ The no-flow outer boundary condition of equation (5) requires that

$\begin{matrix} {A_{0} = {{\frac{C}{4}\sinh\; 2\xi_{e}} > 0.}} & (18) \end{matrix}$ The complete pseudo-steady-state solution for the dimensionless pressure in the reservoir (ξ₁≤ξ≤ξ_(e)) is then

$\begin{matrix} {{p_{D,{PSS}}\left( {\xi,\eta,t_{DL}} \right)} = {B_{0} + {\frac{C}{4}{\xi sinh2\xi}_{e}} - {\sum\limits_{n = 1}^{\infty}{A_{n}\cos\; 2\; n\;{\eta cosh2}\;{n\left( {\xi_{e} - \xi} \right)}}} - {\frac{C}{8}\left( {{\cosh\; 2\xi} + {\cos\; 2\eta}} \right)} - {{Ct}_{DL}.}}} & (19) \end{matrix}$ The constant C is directly related to the fluid production rate. The dimensional flux-density q_(d)(η) on the fracture surface ξ=ξ₁ for the fluid entering the fracture from the reservoir is given by

                                    (20) $\begin{matrix} {{{q_{d}(\eta)} = {\frac{\kappa}{\mu}\frac{p_{i,d}}{L\sqrt{{\sinh^{2}\xi_{1}} + {\sin^{2}\eta}}}\frac{\partial p_{D}}{\partial\xi}}}}_{\xi\; = \;\xi_{1}} \\ {= {\frac{\kappa}{\mu}\frac{p_{i,d}}{L\sqrt{{\sinh^{2}\xi_{1}} + {\sin^{2}\eta}}}}} \\ {\left\lbrack {{\frac{C}{4}\left( {{\sinh\; 2\;\xi_{e}} - {\sinh\; 2\;\xi_{1}}} \right)} + {\sum\limits_{n = 1}^{\infty}{2n\; A_{n}\cos\; 2n\;\eta\;\sinh\; 2{n\left( {\xi_{e} - \xi_{1}} \right)}}}} \right\rbrack.} \end{matrix}$ Therefore the dimensional production-rate for a bi-wing fractured-well is

$\begin{matrix} {{Q_{d} = {{4h{\int_{0}^{\pi/2}{{q_{d}(\eta)}L\sqrt{{\sinh^{2}\xi_{1}} + {\sin^{2}\eta}}d\;\eta}}} = {\frac{\pi\;\kappa\; h\; p_{i,d}}{2\;\mu}{C\left( {{\sinh\; 2\;\xi_{e}} - {\sinh\; 2\;\xi_{1}}} \right)}}}},} & (21) \end{matrix}$ where h is the formation thickness. Thus, the dimensionless parameter C is related to the well production rate by

$\begin{matrix} {C = {\frac{2\;\mu}{\pi\;\kappa\; h\; p_{i,d}}{\frac{Q_{d}}{{\sinh\; 2\;\xi_{e}} - {\sinh\; 2\;\xi_{1}}}.}}} & (22) \end{matrix}$ The Pseudo-Steady State Pressure Profile in the Fracture

The fracture pressure of equation (7) can be written as

$\begin{matrix} {{\frac{\partial^{2}{p_{{fD},{PSS}}\left( {\eta,t_{DL}} \right)}}{\partial\eta^{2}} + {\frac{2}{F_{E}}\left\lbrack {{\frac{C}{4}\left( {{\sinh\; 2\;\xi_{e}} - {\sinh\; 2\;\xi_{1}}} \right)} + {\sum\limits_{n = 1}^{\infty}{2n\; A_{n}\cos\; 2n\;\eta\;\sinh\; 2{n\left( {\xi_{e} - \xi_{1}} \right)}}}} \right\rbrack}} = 0.} & (23) \end{matrix}$ There is no-flow across the x-axis due to symmetry, η=0: ∂p _(fD,PSS)/∂η=0,0≤ξ≤ξ₁.  (24) The pseudo-steady state property also holds for the pressure inside the fracture,

$\begin{matrix} {\frac{\partial{p_{{fD},{PSS}}\left( {\eta,t_{DL}} \right)}}{\partial t_{DL}} = {- {C.}}} & (25) \end{matrix}$ Integration of equation (23) subject to equations (24) and (25) gives the dimensionless fracture pressure

$\begin{matrix} {{{p_{{fD},{PSS}}\left( {\eta,t_{DL}} \right)} = {{- {\frac{2}{F_{E}}\left\lbrack {{\frac{C}{8}\left( {{\sinh\; 2\;\xi_{e}} - {\sinh\; 2\;\xi_{1}}} \right)\eta^{2}} - {\sum\limits_{n = 1}^{\infty}{\frac{A_{n}}{2n}\cos\; 2\; n\;\eta\;\sinh\; 2{n\left( {\xi_{e} - \xi_{1}} \right)}}}} \right\rbrack}} - {Ct}_{DL} + \overset{\sim}{C}}},} & (26) \end{matrix}$ where {tilde over (C)} is an integration constant.

The dimensionless pressure at the well, which is unknown for the PSS solution, is given by

$\begin{matrix} {p_{{wD},{PSS}} = {{p_{{fD},{PSS}}\left( {{\pi\text{/}2},t_{DL}} \right)} = {{- {\frac{2}{F_{E}}\left\lbrack {{\frac{\pi^{2}C}{32}\left( {{\sinh\; 2\;\xi_{e}} - {\sinh\; 2\xi_{1}}} \right)} - {\sum\limits_{n = 1}^{\infty}{\left( {- 1} \right)^{n}\frac{A_{n}}{2n}\sinh\; 2{n\left( {\xi_{e} - \xi_{1}} \right)}}}} \right\rbrack}} - {Ct}_{DL} + {\overset{\sim}{C}.}}}} & (27) \end{matrix}$ Determination of the Coefficients

The coefficients A_(n) and the constant C in the solution for the pressure can be obtained by matching the reservoir pressure on the fracture surface with the fracture pressure and an application of the material balance equation. Because the fracture is narrow and ξ₁ is very small, we set ξ₁≈0 in all calculations below.

Pressure Matching on the Fracture Surface

On the fracture surface, ξ=ξ₁=0, the reservoir pressure and the fracture pressure must match, p _(D,PSS)(0,η,t _(DL))=p _(fD,PSS)(η,t _(DL)). This leads to equation (28):

$\begin{matrix} {{B_{0} - {\sum\limits_{n = 1}^{\infty}{A_{n}\cos\; 2n\;{\eta cosh2}\; n\;\xi_{e}}} - {\frac{C}{8}\left( {1 + {\cos\; 2\eta}} \right)} - {Ct}_{DL}} = {{- {\frac{2}{F_{E}}\left\lbrack {{\frac{C}{8}\eta^{2}\sinh\; 2\;\xi_{e}} - {\sum\limits_{n = 1}^{\infty}{\frac{A_{n}}{2n}\cos\; 2n\;{\eta sinh2}\; n\;\xi_{e}}}} \right\rbrack}} - {Ct}_{DL} + {\overset{\sim}{C}.}}} & (28) \end{matrix}$

It is observed that, without the infinite series in the reservoir pressure, it would not be possible to match the η² term from the fracture pressure. Pressure matching can be accomplished by simply expanding η² as a cosine series.

The Material Balance Equation

The dimensionless reservoir pressure drawdown Δp_(mD) is defined as

$\begin{matrix} {{\Delta\; p_{D}} = {\frac{2{\pi\kappa}\; h}{\mu\; Q_{d}} = {\left\lbrack {p_{i,d} - p_{d}} \right\rbrack = {{\frac{2{\pi\kappa}\; h\; p_{i,d}}{\mu\; Q_{d}}\left\lbrack {1 - p_{D}} \right\rbrack}.}}}} & (29) \end{matrix}$ For a closed reservoir and constant rate production, the material balance equation provides a simple relation between the reservoir average pressure drawdown and time,

$\begin{matrix} {{{\Delta\;{\overset{\_}{p}}_{D}} = {{\frac{2\pi\;\kappa\; h}{\mu\; Q_{d}}\left\lbrack {p_{i,d} - {{\overset{\_}{p}}_{d}(t)}} \right\rbrack} = {2\pi\; t_{DA}}}},} & (30) \end{matrix}$ where p _(d)(t) is the reservoir volume-averaged pressure

$\begin{matrix} {{{{\overset{\_}{p}}_{d}(t)} = {\frac{1}{V}{\int_{V}{p_{d}\ {dv}}}}},} & (31) \end{matrix}$ V being the reservoir volume; and t_(DA) is the dimensionless time defined in terms of the draining area, A=V/h=πL² sin h2ξ_(e)/2,

$\begin{matrix} {t_{DA} = {\frac{\kappa\; t}{\mu\; c_{t}\phi\; A} = {\frac{2}{\pi}{\frac{t_{DL}}{{\pi sinh2\xi}_{e}}.}}}} & (32) \end{matrix}$ Computing the reservoir average pressure using the solution of equation (19) and utilizing the relation between the constant C and the production rate Q_(d) of equation (22), the material balance equation (30) becomes,

$\begin{matrix} {{{\int_{A}{\left\lbrack {1 - B_{0} - \ {\frac{C}{4}{\xi sinh2\xi}_{e}} + {\sum\limits_{n = 1}^{\infty}{A_{n}\cos\; 2n\;{\eta cosh2}\;{n\left( {\xi_{e} - \xi} \right)}}} + {\frac{C}{8}\left( {{\cosh\; 2\xi} + {\cos\; 2\eta}} \right)}} \right\rbrack{dA}}} = 0},\quad} & (33) \end{matrix}$ which leads to

$\begin{matrix} {B_{0} = {1 - {\frac{C}{4}\left( {{\xi_{e}\sinh\; 2\xi_{e}} - \frac{{\cosh\; 2\xi_{e}} - 1}{2}} \right)} + {\frac{C}{16}\cosh\; 2\xi_{e}} - {\frac{A_{1}}{2}.}}} & (34) \end{matrix}$ Matching Fourier coefficients in equation (28) then gives

$\begin{matrix} {{\overset{\sim}{C} = {1 - \frac{C}{4} - {\frac{C}{4}\xi_{e}\sinh\; 2\xi_{e}} + {\frac{3C}{16}\cosh\; 2\xi_{e}} - \frac{A_{1}}{2} + {\frac{C}{F_{E}}\frac{\pi^{2}}{48}\sinh\; 2\xi_{e}}}},} & (35) \\ {\mspace{79mu}{{A_{1} = {- {\frac{C}{{\cosh\; 2\xi_{e}} + \frac{\sinh\; 2\xi_{e}}{F_{E}}}\left\lbrack {\frac{1}{8} + {\frac{1}{F_{E}}\frac{\sinh\; 2\xi_{e}}{4}}} \right\rbrack}}},}} & (36) \\ {\mspace{79mu}{{A_{n} = {\frac{C}{4}\frac{\left( {- 1} \right)^{n}}{n}\frac{\sinh\; 2\xi_{e}}{{\sinh\; 2n\;\xi_{e}} + {n\; F_{E}\cosh\; 2n\;\xi_{e}}}}},{n \geq 2.}}} & (37) \end{matrix}$ Thus, the pressure in the reservoir and the fracture are completely determined. In particular, the dimensionless pressure drawdown in the reservoir is given by

$\begin{matrix} \begin{matrix} {{{\Delta\;{p_{D,{PSS}}\left( {\xi,\eta,t_{DL}} \right)}} = {{2\pi\; t_{DA}} + \xi_{e} - {\xi++}}}} \\ {{\frac{1}{\sinh\; 2\xi_{e}}\begin{bmatrix} {{- \frac{{3\cosh\; 2\xi_{e}} - 2}{4}} + {2a_{1}} +} \\ {{4a_{1}\cos\; 2{{\eta cosh2}\left( {\xi_{e} - \xi} \right)}} +} \\ \frac{{\cosh\; 2\xi} + {\cos\; 2\eta}}{2} \end{bmatrix}}, +} \\ {\sum\limits_{n = 2}^{\infty}{\frac{\left( {- 1} \right)^{n}}{n}\frac{\cos\; 2n\;{\eta cosh}\; 2{n\left( {\xi_{e} - \xi} \right)}}{{\sinh\; 2n\;\xi_{e}} + {n\; F_{E}\cosh\; 2n\;\xi_{e}}}}} \end{matrix} & (38) \\ {\mspace{79mu}{where}} & \; \\ {\mspace{79mu}{a_{1} = {\frac{A_{1}}{C} = {{- \frac{1}{8}}{{\frac{1}{{\cosh\; 2\xi_{e}} + \frac{\sinh\; 2\xi_{e}}{F_{E}}}\left\lbrack {1 + {\frac{2}{F_{E}}\sinh\; 2\xi_{e}}} \right\rbrack}.}}}}} & (39) \end{matrix}$ Shank's transformation can be used to accelerate the convergence of the infinite series in equation (38).

The dimensionless pressure drawdown at the well is given by

$\begin{matrix} {{\Delta\; p_{{wD},{PSS}}} = {{2\;\pi\; t_{DA}} + \xi_{e} + \frac{1}{\sinh\mspace{11mu} 2\;\xi_{e}} - {\frac{3}{4}\coth\mspace{11mu} 2\;\xi_{e}} + \frac{2\; a_{1}}{\sinh\mspace{11mu} 2\xi_{e}} + {{\frac{1}{F_{E}}\left\lbrack {\frac{\pi^{2}}{6} + {4a_{1}} - {\sum\limits_{n = 2}^{\infty}\;{\frac{1}{n^{2}}\frac{1}{1 + {n\; F_{E}\coth\mspace{11mu} 2\; n\;\xi_{e}}}}}} \right\rbrack}.}}} & (40) \end{matrix}$ Thus, an explicit expression for the pseudo-steady state constant b_(D,PSS) is given by

$\begin{matrix} {b_{D,{PSS}} = {\xi_{e} + \frac{1}{\sinh\mspace{11mu} 2\;\xi_{e}} - {\frac{3}{4}\coth\mspace{11mu} 2\;\xi_{e}} + \frac{2\; a_{1}}{\sinh\mspace{11mu} 2\;\xi_{e}} + {{\frac{1}{F_{E}}\left\lbrack {\frac{\pi^{2}}{6} + {4a_{1}} - {\sum\limits_{n = 2}^{\infty}\;{\frac{1}{n^{2}}\frac{1}{1 + {n\; F_{E}\coth\mspace{11mu} 2\; n\;\xi_{e}}}}}} \right\rbrack}.}}} & (41) \end{matrix}$ In addition, the productivity index (PI) and the dimensionless productivity index (J_(D)) for the pseudo-steady state flow is given by

$\begin{matrix} {\begin{matrix} {J_{PSS} = {\frac{Q_{d}}{{\overset{\_}{P}}_{d} - P_{w,d}} = {\frac{\kappa\; h}{\mu}\frac{2\;\pi}{b_{D,{PSS}}}}}} \\ {J_{D,{PSS}} = {{\frac{\mu}{2\;{\pi\kappa}\; h}J_{PSS}} = \frac{1}{b_{D,{PSS}}}}} \end{matrix}.} & (42) \end{matrix}$ The dimensionless productivity index J_(D), or the effective wellbore radius, can be used to characterize the productivity of unfractured and fractured wells. For example, J_(D,PSS) can be used for fracture design.

FIG. 3 is a flow chart illustrating an example method for computing an analytical solution for pseudo-steady state flow for a vertically fractured well with finite fracture conductivity in a closed reservoir according to some embodiments of the disclosure. A method 300 may begin at block 302 with receiving one or more shape factors corresponding to a geometrical shape of a hydraulically fractured well reservoir. The data received at block 302 may be received through, for example, an input device or local storage coupled to a processor or may be received through a network communication from a remote data store or remote input device. Examples of the one or more shape factors include ellipse focal distance/fracture half-length, formation thickness, dimensionless elliptical fracture conductivity, wellbore radius, radius of circular drainage boundary, reservoir volume, fracture width at the wellbore, elliptical coordinates, elliptical fracture shape, and elliptical reservoir shape.

Then, at block 304, a pseudo-steady state constant may be determined by the processor, such as using equation (41) for the reservoir based on the plurality of shape factors. The determination at block 304 may be performed using one or more elementary functions to obtain an analytical solution and/or without solving Mathieu functions, which can significantly improve the computational speed of the determination in comparison to prior art numerical simulations. Block 304 may alternatively or additionally include a computation of reservoir pressure drawdown from, for example, equation (38).

Next, at block 306, one or more performance parameters of the reservoir may be determined by the processor when the reservoir is operated in a pseudo-steady state with a finite fracture conductivity based on the determined pseudo-steady state constant. Although block 306 describes finite fracture conductivity, infinite fracture conductivity may alternatively be used for determining the performance parameter. Examples of performance parameters include a production decline rate for a reservoir, a total hydrocarbon reserves for a reservoir, an economically-recoverable reserves for a reservoir, the productivity index (PI), and the dimensionless productivity index (J_(D)). Using the performance parameters, decisions as to explore and produce from certain reservoirs may be made, and the improved information available from the pseudo-steady state analysis of the method 300 may increase profitability of the production from selected reservoirs. The one or more performance parameters or the pseudo-steady state constant may be stored in local or remote storage, output to a display screen, or communicated to another device through a network communications connection. Additional computations or decisions may be performed using the performance parameter, such as decisions relating to the production of hydrocarbons from a particular reservoir.

The specific features of the method 300 for determining a pseudo-steady state constant and a performance parameter from that constant results in a specific process for evaluating reservoirs using particular information and techniques. Analysis of reservoirs using the method 300 results in a technological improvement over the prior art numerical solutions, which are tedious simulations to process. The method 300 thus describes a process specifically designed to achieve an improved technological result of decreased computational time and increased computational accuracy in the conventional industry practice of determining performance from reservoirs. Furthermore, the method 300, and particularly block 304, describes a new analytical solution for calculation of parameters related to a reservoir that differs from conventional industry solutions.

FIG. 4 are graphs of a pseudo-steady state constant computation as a function of ξ_(e) calculated according to some embodiments of the disclosure. The pseudo-steady state constant b_(D,PSS)(ξ_(e),F_(E)) is plotted against ξ_(e) for F_(E)=1, 2, 5, 10, 20, 1000 in lines 402, 404, 406, 408, 410, and 412 of FIG. 4, respectively. It is observed that for large ξ_(e), the slope ∂b_(D,PSS)/∂ξ_(e) becomes one, regardless of the value of the fracture conductivity F_(E).

Comparison with Existing Results

The analytical solution obtained in the present work is exact and general under the assumptions adopted, and the solution is valid for both infinite and finite fracture conductivities. A comparison between this new analytical solution and presently known results is provided below.

Infinite Fracture Conductivity.

One conventional pseudo-steady state solution for the case of infinite fracture conductivity, F_(E)→∞ shows that the dimensionless pressure drawdown in the reservoir:

$\begin{matrix} {{{\Delta\; p_{D,{Prats}}} = {{2\; t_{D,{Prat}}} + \xi_{e} + \frac{1}{2\mspace{11mu}\sinh\mspace{11mu} 2\;\xi_{e}} - {\frac{3}{4}\coth\mspace{11mu} 2\;\xi_{e}} - \frac{1}{2\mspace{11mu}\sinh\mspace{11mu} 4\xi_{e}} - \xi - {\frac{1}{\sinh\mspace{11mu} 4\xi_{e}}\cosh\mspace{11mu} 2\left( {\xi_{e} - \xi} \right)\cos\mspace{11mu} 2\;\eta} + \frac{{\cosh\mspace{11mu} 2\xi} + {\cos\mspace{11mu} 2\eta}}{2\mspace{11mu}\sinh\mspace{11mu} 2\xi_{e}}}},} & (43) \end{matrix}$ where dimensionless time t_(D,Prat) is defined as related to t_(DA) by (after a correction to a missing factor ϕ in their definition):

$\begin{matrix} {t_{D,{Prats}} = {\frac{\pi\;\kappa\; t}{{\mu\phi}\; c_{t}A} = {\pi\;{t_{DA}.}}}} & (44) \end{matrix}$ Thus, the dimensionless reservoir pressure drawdown from is:

$\begin{matrix} {{\Delta\; p_{D,{Prats}}} = {{2\pi\; t_{DA}} + \xi_{e} + \frac{1}{2\mspace{11mu}\sinh\mspace{11mu} 2\;\xi_{e}} - {\frac{3}{4}\coth\mspace{11mu} 2\;\xi_{e}} - \frac{1}{2\mspace{11mu}\sinh\mspace{11mu} 4\xi_{e}} - \xi - \frac{\cosh\mspace{11mu} 2\left( {\xi_{e} - \xi} \right)\cos\mspace{11mu} 2\;\eta}{\sinh\mspace{11mu} 4\xi_{e}} + {\frac{{\cosh\mspace{11mu} 2\xi} + {\cos\mspace{11mu} 2\eta}}{2\mspace{11mu}\sinh\mspace{11mu} 2\xi_{e}}.}}} & (45) \end{matrix}$

From equation (38), for infinite fracture conductivity, F_(E)→∞, the infinite sum becomes zero, and

$a_{1} = {{- \frac{1}{8}}{\frac{1}{\cosh\mspace{11mu} 2\xi_{e}}.}}$ Thus, the reservoir dimensionless pressure drawdown becomes

${{\Delta\;{p_{D,{PSS}}\left( {\xi,\eta,t_{DL}} \right)}} = {{2\pi\; t_{DA}} + \xi_{e} - {{\xi++}{\frac{1}{\sinh\mspace{11mu} 2\xi_{e}}\begin{bmatrix} {{- \frac{{3\mspace{11mu}\cosh\mspace{11mu} 2\xi_{e}} - 2}{4}} - {\frac{1}{4}\frac{1}{\cosh\mspace{11mu} 2\xi_{e}}}} \\ {\;{{{- \frac{1}{2}}\frac{1}{\cosh\mspace{11mu} 2\xi_{e}}\cos\mspace{11mu} 2\eta\mspace{11mu}\cosh\mspace{11mu} 2\left( {\xi_{e} - \xi} \right)} + \frac{{\cosh\mspace{11mu} 2\xi} + {\cos\mspace{11mu} 2\eta}}{2}}} \end{bmatrix}}}}},$ which is identical to the prior art result of equation (45).

Similarly, the pressure drawdown at the well from our solution becomes

$\begin{matrix} {{{\Delta\; p_{D,{PSS}}2\pi\; t_{DA}} + \xi_{e} - \frac{\left( {{3\mspace{11mu}\cosh\mspace{11mu} 2\xi} - 1} \right)\left( {{\cosh\mspace{11mu} 2\xi} - 1} \right)}{4\mspace{11mu}\sinh\mspace{11mu} 2\xi_{e}\mspace{11mu}\cosh\mspace{11mu} 2\xi_{e}}},} & (46) \end{matrix}$ which gives the pseudo-steady state constant for the case of infinite fracture conductivity

$\begin{matrix} {b_{D,{PSS}} = {\xi_{e} - {\frac{\left( {{3\mspace{11mu}\cosh\mspace{11mu} 2\xi_{e}} - 1} \right)\left( {{\cosh\mspace{11mu} 2\xi_{e}} - 1} \right)}{4\mspace{11mu}\sinh\mspace{11mu} 2\xi_{e}\mspace{11mu}\cosh\mspace{11mu} 2\xi_{e}}.}}} & (47) \end{matrix}$

In summary, in the limit of infinite fracture conductivity, an analytical solution according to embodiments described herein matches a conventional solution for infinite fracture conductivity. This demonstrates that the analytical model is correct, and that at least one specific calculation matches a result from a conventional model.

Finite Fracture Conductivity.

For finite fracture conductivity, the pseudo-steady state constant also depends on the dimensionless fracture conductivity F_(E):b_(D,PSS)=b_(D,PSS)(ξ_(e),F_(E)). b_(D,PSS)(ξ_(e),F_(E)) has been computed in the prior art for selected sets of ξ_(e),F_(E) by subtracting 2πt_(DA) from numerical simulation results for large times. This procedure involves numerical manipulation of the Mathieu functions in the Laplace transform space as well as numerical inversion; and it is tedious and time-consuming, as noted by these authors. A nonlinear-regression may be applied to fit such numerical results into an empirical formula for b_(D,PSS)(ξ_(e),F_(E))

$\begin{matrix} {{{b_{D,{PSS}}\left( {\xi_{e},F_{E}} \right)} = {{1.00146\;\xi_{e}} + {0.0794849\; e^{- \xi_{e}}} - {0.16703\; u} + \frac{A}{B} - 0.754772}},} & (48) \end{matrix}$ with u=ln F_(E), A=a ₁ +a ₂ u+a ₃ u ² +a ₄ u ³ +a ₅ u ⁴ ,B=b ₁ +b ₂ u+b ₃ u ² +b ₄ u ³ +b ₅ u ⁴, a ₁=−4.7468,b ₁=−2.4941, a ₂=36.2492,b ₂=21.6755, a ₃=55.0998,b ₃=41.0303, a ₄=−3.98311,b ₄=−10.4793, a ₅=6.07102,b ₅=5.6108.

However, there are some apparent inconsistency and problematic issues with equation (48): (i) the formula cannot re-produce certain tabulated results of the prior art; (ii) equation (48) can give rise to negative values of b_(D,PSS) when F_(E) becomes large; and it does not converge to the exact result of the prior art for infinite fracture conductivity; (iii) when the empirical equation (48) is compared to the disclosed analytical solution for b_(D,PSS)(ξ_(e),F_(E)) in equation (41), it is immediately obvious that the coefficient for the linear term ξ_(e) in equation (48) must be “1.0”, instead of “1.00146.”

The results of the analytical solution of equation (41) are computed and compared to corresponding values b_(D,PSS)(ξ_(e),F_(E)) for the parameter sets as known in the prior art. The results are shown below in our Table 1, where the results of prior art are listed in the parentheses for comparison. It is seen that the numerically computed values from the prior art generally agree very well with the described analytical solution.

TABLE 1 Values of b_(D,PSS)(ξ_(e),F_(E)) from the analytical solution. Values in the parentheses are those of prior art numerical simulations. ξ_(e) F_(E) = 1 F_(E) = 10 F_(E) = 100 F_(E) = 1000 0.25 0.849411 (0.8481)  0.213087 (0.2150) 0.130127 (0.1306) 0.121565 (0.1220) 0.50 0.989853 (0.9902)  0.333336 (0.3337) 0.239246 (0.2396) 0.229383 (0.2298) 0.75 1.16694 (1.1671) 0.460557 (0.4609) 0.353713 (0.3540) 0.342402 (0.3426) 1.00  1.3632 (1.3627) 0.610541 (0.6109) 0.493289 (0.4936) 0.480809 (0.4812) 1.25 1.57305 (1.5733) 0.787704 (0.7880) 0.663153 (0.6634) 0.649857 (0.6501) 1.50 1.79635 (1.7963) 0.988308 (0.9884) 0.858987 (0.8591) 0.845162 (0.8453) 1.75 2.02893 (2.0293)  1.20624 (1.2067)  1.07391 (1.0743)  1.05975 (1.0602) 2.00 2.26787 (2.2682)  1.43597 (1.4363)  1.30178 (1.3021)  1.28741 (1.2877) 3.00 3.25252 (3.2529)  2.40795 (2.4084)  2.27122 (2.2716)  2.25658 (2.2570) 4.00 4.25038 (4.2503)  3.40407 (3.4040)  3.26699 (3.2669)  3.25231 (3.2522) 5.00 5.25009 (5.2486)  4.40354 (4.4021)  4.26642 (4.2649)  4.25173 (4.2502) Discussions

Analytical solutions for the reservoir pressure drawdown Δp_(D,PSS) and the pseudo-steady state constant b_(D,PSS) are given by equations (38) and (41), respectively. These expressions may be exact for fully-penetrating hydraulically fractured vertical wells producing from a closed reservoir approximated as having an elliptical shape, and the solutions are valid for both finite and infinite fracture conductivities. As a result of the analytical solution described in embodiments of the disclosure herein, tedious and time consuming numerical simulations for obtaining pseudo-steady state solutions for fractured wells are no longer necessary for such reservoirs. For a fractured-well in a large circular reservoir with a radius R, these formulas can be readily applied with ξ_(e)=ln(2R/L), because a large circle and an ellipse with large ξ_(e) are nearly identical. It is also possible to extend the expression for the pseudo-steady state constant b_(D,PSS) to a fractured-well in reservoirs of different geometrical shapes using an equivalent elliptical parameter ξ_(e) based on the reservoir drainage area or shape factors.

Conclusions

An exact analytical solution for pseudo-steady state productive flow from a fully-penetrating hydraulically fractured vertical well with finite fracture conductivity in a closed reservoir approximated as having an elliptical shape is rigorously derived. The solution agrees with prior art solutions in the limit of infinite fracture conductivity, and it agrees with the numerical results of the prior art for finite fracture conductivity. The analytical solution is exact, general and expressed in terms of elementary functions; it is simple and easy to evaluate; and it completely eliminates the need of performing numerical simulation for obtaining pseudo-steady state solution for a vertically fractured well in such a reservoir. Simple expressions for the pseudo-steady state constant and the dimensionless productivity index are described above. The solution may also be used to generate approximate analytical solutions for pseudo-steady state flow from a fractured-well in reservoirs of different geometrical shapes.

Advantages of Embodiments of the Invention

An exact analytical solution in the physical variable space for pseudo-steady state production from a vertically fractured well with finite fracture conductivity in an elliptical reservoir is obtained from this work. The solution is expressed in terms of elementary functions and it yields a simple, exact expression for the pseudo-steady state constant b_(D,PSS) and the dimensionless productivity index J_(D,PSS). This is the first time that an exact analytical solution has been obtained for pseudo-steady state flow for a fractured well with finite conductivity.

Some advantages resulting from this analytical solution are listed below:

(1) It eliminates the need to perform time-consuming numerical simulation in order to obtain the pseudo-steady state constant b_(D,PSS) b_(D,PSS) and the dimensionless productivity index J_(D,PSS) for fractured wells in elliptical reservoirs, and it shortens the required computing time from hours/days to seconds;

(2) By introducing suitable shape factors, the solution can be used to obtain approximate expressions for the pseudo-steady state constant b_(D,PSS) for fractured wells in reservoirs of other geometrical shapes;

(3) The solution can be readily adopted for use with production decline models and simulators for estimating total hydrocarbon reserves in-place as well as economically recoverable reserves;

(4) The solution can be used for optimal fracture design so that the production is optimized;

(5) The solution can be used as a benchmark to measure the accuracy of various numerical simulators; and

(6) The techniques used in certain embodiments of the disclosure (such as hyperbolic functions and Fourier series expansions) circumvent the cumbersome Mathieu functions commonly used in studying production from fractured wells, and these techniques can be adopted for much wider use in studying similar problems.

Important of Pseudo Steady-State Flow Analysis

The duration of the pseudo-steady state flow and its productive performance largely determines the cumulative production of hydrocarbon form a well. The duration of pseudo-steady state flow is determined by how fast the bottom-hole flowing pressure decreases to the lowest permissible well pressure (critical pressure). Thus it is paramount to know the change of the wellbore pressure with time. The pressure drawdown (pressure drop from the initial reservoir pressure) at the wellbore is commonly expressed in dimensionless form as Δp _(wD,PSS)=2πt _(D4) +b _(D,PSS),  (49) where t_(DA) is the drainage area based dimensionless time, and b_(D,PSS) is the so-called pseudo-steady state constant which depends on the reservoir model as well as the well/reservoir configuration. Thus, two parameters determine the duration of the pseudo-steady state flow period: the time-rate of decline, which is determined by the production rate, and the pseudo-steady state constant.

The productive performance of a well is measured by the productivity index, J, which is the amount of hydrocarbon produced per unit drop in the reservoir average pressure. For pseudo-steady state flow, the productivity index is inversely proportional to the pseudo-steady state constant b_(D,PSS),

$\begin{matrix} {{J = {\frac{\kappa\; h}{\mu}\frac{2\;\pi}{b_{D,{PSS}}}}},} & (50) \end{matrix}$ where κ,μ,h are the reservoir permeability, hydrocarbon viscosity, and hydrocarbon bearing formation thickness, respectively. Thus, the productivity of a well during the pseudo-steady state period is completely determined by the pseudo-steady state constant b_(D,PSS).

Furthermore, pseudo-steady state solution has been often used in the production rate decline analysis because pseudo-steady state is the flow regime immediate preceding the production rate decline period (as shown in FIG. 1). Production rate decline can be used for estimating the hydrocarbon reserves in place and for assessing the economically recoverable amount of hydrocarbon from a reservoir.

The pseudo-steady state flow analysis can be used to improve production from reservoirs, because: Pseudo-steady state flow can impact the cumulative production of hydrocarbon from a well; the productive performance of a well can be assessed by evaluating the productivity of the well during the pseudo-steady state flow, which is determined by the value of the pseudo-steady state constant b_(D,PSS); Pseudo-steady state flow can be used for estimating the total reserves in place in a reservoir; and Pseudo-steady state flow can be used for estimating the economically recoverable amount of hydrocarbon from a reservoir.

Implementation

Computations described in the embodiments above may be executed on any suitable processor-based device including, without limitation, personal data assistants (PDAs), tablet computers, smartphones, computer game consoles, and multi-processor servers. Moreover, the systems and methods of the present disclosure may be implemented on application specific integrated circuits (ASIC), very large scale integrated (VLSI) circuits, or other circuitry.

If implemented in firmware and/or software, the functions described above may be stored as one or more instructions or code on a computer-readable medium. Examples include non-transitory computer-readable media encoded with a data structure and computer-readable media encoded with a computer program. Computer-readable media includes physical computer storage media. A storage medium may be any available medium that can be accessed by a computer. By way of example, and not limitation, such computer-readable media can comprise RAM, ROM, EEPROM, CD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other medium that can be used to store desired program code in the form of instructions or data structures and that can be accessed by a computer. Disk and disc includes compact discs (CD), laser discs, optical discs, digital versatile discs (DVD), floppy disks and blu-ray discs. Generally, disks reproduce data magnetically, and discs reproduce data optically. Combinations of the above should also be included within the scope of computer-readable media.

In addition to storage on computer readable medium, instructions and/or data may be provided as signals on transmission media included in a communication apparatus. For example, a communication apparatus may include a transceiver having signals indicative of instructions and data. The instructions and data are configured to cause one or more processors to implement the functions outlined in the claims.

Although the present disclosure and its advantages have been described in detail, it should be understood that various changes, substitutions and alterations can be made herein without departing from the spirit and scope of the disclosure as defined by the appended claims. Moreover, the scope of the present application is not intended to be limited to the particular embodiments of the process, machine, manufacture, composition of matter, means, methods and steps described in the specification. As one of ordinary skill in the art will readily appreciate from the present invention, disclosure, machines, manufacture, compositions of matter, means, methods, or steps, presently existing or later to be developed that perform substantially the same function or achieve substantially the same result as the corresponding embodiments described herein may be utilized according to the present disclosure. Accordingly, the appended claims are intended to include within their scope such processes, machines, manufacture, compositions of matter, means, methods, or steps. 

What is claimed is:
 1. A method for increasing production of hydrocarbon from a hydraulically fractured well reservoir, the method comprising: receiving over an electronic network, from a remote input device and at a processor configured for determining parameters of hydrocarbon production for the hydraulically fractured well reservoir, a plurality of shape factors corresponding to a geometrical shape of a hydraulically fractured well reservoir; determining, by the processor, a pseudo-steady state constant for the reservoir based, at least in part, on an analytical solution involving the plurality of shape factors; determining, by the processor, a performance parameter of the reservoir when operated in a pseudo-steady state based on the determined pseudo-steady state constant; and selecting, based on the performance parameter, the reservoir for increasing production of hydrocarbon, wherein the pseudo-steady state constant is determined, by the processor, according to the following equation: $b_{D,{PSS}} = {\xi_{e} + \frac{1}{\sinh\mspace{11mu} 2\;\xi_{e}} - {\frac{3}{4}\coth\mspace{11mu} 2\;\xi_{e}} + \frac{2\; a_{1}}{\sinh\mspace{11mu} 2\;\xi_{e}} + {{\frac{1}{F_{E}}\left\lbrack {\frac{\pi^{2}}{6} + {4a_{1}} - {\sum\limits_{n = 2}^{\infty}\;{\frac{1}{n^{2}}\frac{1}{1 + {n\; F_{E}\coth\mspace{11mu} 2\; n\;\xi_{e}}}}}} \right\rbrack}.}}$
 2. The method of claim 1, further comprising ceasing production of hydrocarbons from the reservoir based on the performance parameter.
 3. The method of claim 1, wherein the step of determining the pseudo-steady state constant is performed without solving Mathieu functions.
 4. The method of claim 1, wherein the step of determining the performance parameter comprises determining total hydrocarbon reserves for the reservoir.
 5. The method of claim 1, wherein the shape factors comprise one or more of ellipse focal distance, fracture half-length, formation thickness, dimensionless elliptical fracture conductivity, wellbore radius, radius of circular drainage boundary, reservoir volume, fracture width at the wellbore, elliptical coordinates, elliptical fracture shape, or elliptical reservoir shape.
 6. The method of claim 1, wherein the determining the performance parameter comprises determining a production decline rate for the reservoir.
 7. The method of claim 1, wherein the determining the performance parameter comprises determining economically recoverable reserves for the reservoir.
 8. The method of claim 1, wherein a wellhead is coupled to a hydraulically fractured vertical well that fully penetrates the reservoir.
 9. A method for increasing a increasing production of hydrocarbon from a hydraulically fractured well reservoir, the method comprising: receiving over an electronic network, from a remote input device and at a processor configured for determining parameters of hydrocarbon production for the hydraulically fractured well reservoir, a plurality of shape factors corresponding to a geometrical shape of a hydraulically fractured well reservoir; determining, by the processor, a pseudo-steady state constant for the reservoir based, at least in part, on an analytical solution involving the plurality of shape factors; determining, by the processor, a performance parameter of the reservoir when operated in a pseudo-steady state based on the determined pseudo-steady state constant; and assessing, based on the performance parameter, an economically recoverable amount of hydrocarbon from the reservoir, wherein the pseudo-steady state constant is determined, by the processor, according to the following equation: $b_{D,{PSS}} = {\xi_{e} + \frac{1}{\sinh\mspace{11mu} 2\;\xi_{e}} - {\frac{3}{4}\coth\mspace{11mu} 2\;\xi_{e}} + \frac{2\; a_{1}}{\sinh\mspace{11mu} 2\;\xi_{e}} + {{\frac{1}{F_{E}}\left\lbrack {\frac{\pi^{2}}{6} + {4a_{1}} - {\sum\limits_{n = 2}^{\infty}\;{\frac{1}{n^{2}}\frac{1}{1 + {n\; F_{E}\coth\mspace{11mu} 2\; n\;\xi_{e}}}}}} \right\rbrack}.}}$
 10. A method for increasing production of hydrocarbon from a hydraulically fractured well reservoir, the method comprising: receiving over an electronic network, from a remote input device and at a processor configured for determining parameters of hydrocarbon production for the hydraulically fractured well reservoir, a plurality of shape factors corresponding to a geometrical shape of a hydraulically fractured well reservoir; determining, by the processor, a pseudo-steady state constant for the reservoir based, at least in part, on an analytical solution involving the plurality of shape factors; determining, by the processor, a performance parameter of the reservoir when operated in a pseudo-steady state based on the determined pseudo-steady state constant; and estimating, by the processor and based on the performance parameter, the hydrocarbon reserves of the reservoir, wherein the pseudo-steady state constant is determined, by the processor, according to the following equation: $b_{D,{PSS}} = {\xi_{e} + \frac{1}{\sinh\mspace{11mu} 2\;\xi_{e}} - {\frac{3}{4}\coth\mspace{11mu} 2\;\xi_{e}} + \frac{2\; a_{1}}{\sinh\mspace{11mu} 2\;\xi_{e}} + {{\frac{1}{F_{E}}\left\lbrack {\frac{\pi^{2}}{6} + {4a_{1}} - {\sum\limits_{n = 2}^{\infty}\;{\frac{1}{n^{2}}\frac{1}{1 + {n\; F_{E}\coth\mspace{11mu} 2\; n\;\xi_{e}}}}}} \right\rbrack}.}}$ 